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Qubit manipulation is fundamental to quantum computing. However, the stray coupling be- 
tween the qubit and the environment usually present in real experimental systems hinders the 
high-precision control necessary for a quantum computer with error rates below the quantum error 
correction threshold. We consider this problem in the context of the singlet-triplet spin qubit in 
a semiconductor double quantum dot system, where the experimenter has precise control over the 
exchange interaction only, while statistical fluctuations in the magnetic field gradient produce drift 
error about the a:-axis of the Bloch sphere. The exchange interaction is restricted to be positive 
and finite. Taking these constraints into account, we design experimentally viable electrical pulse 
sequences that, for small magnetic field gradients, carry out rotations around z while canceling gate 
errors up to the sixth order in the gradient fluctuations, and for large gradients, perform arbitrary 
rotations while canceling the leading order error. 



A quantum computer would permit exponentially 
faster algorithms than an ordinary computer for certain 
important types of problems^. Universal quantum com- 
putation requires the ability to perform an entangling 
two-qubit gate and precise single-qubit rotations around 
two different axes of the Bloch sphere. Semiconductor 
spin qubits potentially have marked advantages in fast 
two-qubit gating and scalability, but suffer an embar- 
rassing difficulty in performing fast single-qubit rotations 
since the strong, high-frequency magnetic fields naively 
required^ heat the sample and are hard to confine to a 
single qubilP. This problem is circumvented by encoding 
the qubit in the low-lying singlet-triplet subspace of a 
two-electron double quantum dot^^. Fast electrical con- 
trol of the exchange coupling, J, via the tilt of the effec- 
tive double-well potential allows sub-ns rotations about 
the z-axis of the Bloch sphere. 

In order to perform arbitrary rotations of such a qubit, 
though, one must introduce a difference, AB, between 
the local magnetic fields at each dot, resulting in rotation 
about the a;-axis of the Bloch sphere. This can be done 
either by pumping a nuclear spin polarization gradienlP 
or by depositing a micromagnet nearby^. However, a 
problem remains: The local magnetic field typically fluc- 
tuates slowly due to second-order nuclear spin flip-flops 
mediated by the hyperfine coupling to the electron spin® ^ 
and due to charge-noise-induced shifts of the double dot 
position in the inhomogeneous field. The resulting un- 
certainty in AB introduces a quasi-static random com- 
ponent to the rotation about the x-axis which, in the 
course of ensemble averaging, leads to rapid decoherence 
of the qubit on the free induction decay timescale of Tj* • 
In the case of quantum memory, dynamical decoupling 
techniques^ ^'^ can be employed to preserve qubit infor- 
mation long beyond Tj*, up to a time T2 at which this 
information is ultimately lost due to dynamical fluctu- 
ations. However, such echo techniques cannot be per- 
formed simultaneously with arbitrary single-qubit rota- 
tions, so gate errors are still dominated by statistical fluc- 



tuations and depend on the ratio of the gate time to T| 
rather than to T2. It is of utmost importance to address 
this problem since fault-tolerant quantum computation 
requires extremely precise single-qubit rotations. Thus 
the task is to find dynamically corrected gate^r^^ ap- 
plicable to singlet-triplet qubits. In this work, we show 
how to perform dynamically corrected single-qubit gates, 
dramatically reducing errors for experimentally relevant 
scenarios (see Figs, [s] and |4]). 

In this Letter, we report theoretical progress in dynam- 
ically correcting errors stemming from fluctuations in the 
magnetic field gradient in the singlet-triplet spin qubit 
system. These errors are not amenable to correction via 
the numerous quantum control techniques developed in 
other fields, such as NMR, since in the present case the 
experimenter has precise control only over the rotation 
rate about the z-axis of the Bloch sphere, and this rate is 
furthermore restricted to be positive and bounded. De- 
spite these strong constraints, we construct simple elec- 
trical pulse sequences that, for small gradients, carry out 
rotations about the z-axis while canceling errors up to 
the sixth order in the gradient fluctuations, and for large 
gradients, carry out arbitrary rotations while canceling 
the leading order error. This represents an important 
step forward in the development of singlet-triplet qubits 
as viable resources for quantum computing. 

We consider the Hamiltonian governing the singlet- 
triplet qubit. 



h Jit) 
H{t) = -o^ + -^o. 



(1) 



with constraints on the parameters impos ed to accou nt 
for the physical realities of the experiment j^ * ^ * ^^ * ^'^ * ^^* ^. 
Here h = qixbAB and fluctuations in AS are much 
slower than typical gate times so that h = hQ+6h with Hq 
a known constant and Sh a random, unknown constant. 
(T2* is inversely proportional to the width of the statis- 
tical distribution from which 6h is drawn.) The qubit is 
manipulated via the electrically controlled exchange cou- 
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FIG. 1: (Color online.) Example of the five-piece SUPCODE 
pulse canceling gate error up to the fourth order in 5h for 
ho = 0. The inset shows a schematic profile of such a pulse as 
a function of time. The time duration of each piece is denoted 
by ti, t2, and ts. In the main panel we show the values of 
ti (red solid line), t2 (blue dashed line), and (black dotted 
line) as functions of the rotation angle (j). 
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FIG. 2: (Color online.) Example of evolution of a given state 
on the Bloch sphere under the five-piece rotation by 7r/2 about 
z. Here ho = and for the purpose of illustration we have 
taken 5/i/Jmax = 0.05. 



pling, J{t), which is constrained to be positive (except 
in very high magnetic fields, when it is always negative) 
and is restricted in magnitude by the fact that as it is 
increased by tilting the double dot potential, the qubit 
information starts to transfer from the spin state into 
the charge state of the double dot and is hence subject 
to rapid decoherence due to electrical noise, so that one 
must impose a limit, < J(t) < J,nax- The constraint to 
positive rotations is also relevant to exchange-only coded 
qubitPHSl. 

Our goal is to design a pulse in J{t) that respects these 
constraints and performs a given rotation in a way that is 
insensitive to Sh. We dub the resulting composite pulse 
sequence SUPCODE: Soft Uniaxial Positive Control for 
Orthogonal Drift Error. We emphasize that although the 
singlet-triplet qubit is one of the most experimentally ad- 
vanced paths towards a scalable quantum computer, the 
restricted control available does not permit application of 
some elegant methods of quantum control developed re- 
cently. In particular, the positivity of J(t) precludes the 
prescriptions of Refs. I14|16[ and while later related work 
incorporates a positivity constraint^, it is nonetheless 
precluded by the lack of real-time control over h. Below 
we consider two cases: /iq = 0, which is directly rel- 
evant to experiments with unpumped nuclear spins and 
no micromagnelPI^ where only rotations about the 2-axis 
are desired {e.g., for spin echo), and ho ~ Jmax, which 
is directly relevant to experiments with pumped nuclear 
spin^ or a nearby micromagnelP where full single-qubit 
control is desired. 

For ho = 0, we assume for simplicity that the pulse is of 
a binary form where J{t) alternates between its extremal 
values of and Jmax- Expanding the evolution opera- 



tor of the system, U [T) — T|exp —i dtH (t) |, in 

powers of Sh/J^i^^, we choose the time duration of each 
segment such that the zeroth order term is a rotation 
about z by the desired angle and one or more succes- 
sively higher order terms vanish. Details are given in 
the Supplementary Information, where we also rigorously 
show that there is no pulse which cancels all higher or- 
ders for finite Jmax- Thus, we find three-, five-, seven-, 
and nine-piece SUPCODE pulses that perform arbitrary 
rotations about z while canceling undesired terms up to 
the first, second, second, and third order in 5h/ Jmax, re- 
spectively. An example of the five-piece SUPCODE pulse, 
which cancels first and second order terms in i5/i/Jniax 
in the evolution operator (corresponding to second and 
fourth order terms in the gate error) , is shown in Fig. [T] 
Figure [2] traces the evolution of a particular initial state 
on the Bloch sphere under a five-piece SUPCODE pulse 
designed to perform a 7r/2 rotation about z. 
We define the average error per gate. A, as 



(2) 



where V is the desired operation, U{Tf) is the actual 
evolution operator under the composite pulse of dura- 
tion Tf, and the overlap is averaged over initial states 
{ipi) distributed uniformly over the Bloch sphere. Figure 
|3] shows the average error per gate for the naive one- 
piece pulse and the SUPCODE pulses introduced above. 
We see that for 5/i/Jmax < 10% the SUPCODE pulses 
have markedly less error. This range is relevant to recent 
landmark experiments with GaAa^^ and silicon-basecP^ 
double quantum dots. In the former 6h ~ 8neV and in 
the latter Sh ^ 3neV while for both Jmax could be sev- 
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FIG. 3: (Color online.) Average error per gate for SUPCODE 
pulses at /lo = 0. These figures show the average error per 
gate (defined in the main text) as functions of 5/i/Jmax for 
the rotation around the z-axis by (a) <^o = 7r/2 and (b) 
<j!>o — 1.7tv. The "l-piece" curve is for the uncorrected pulse 
with quadratic leading order errors in 5h. The "3-piece," "5- 
piece," "7-piece," and "9-piece" curves are for the three-, five-, 
seven-, and nine-piece pulse sequences that cancels gate error 
up through second, fourth, fourth, and sixth order in 5h, re- 
spectively. For details of the pulses, see the Supplementary 
Information. 





FIG. 4: (Color online.) Average error per gate for SUPCODE 
pulse at ho/Jn-ia.^ ~ 0.25. This set of figures shows the aver- 
age error per gate (defined in the main text) as functions of 
5/i/Jmax, for rotations by n/2 about the (a) x-axis and (b) z- 
axis. The corrected pulses cancel the leading order error. 

operation /(f; (p) designed such that the error in its im- 
plementation exactly cancels the leading order error in 
the original rotation. We represent the uncorrected ro- 
tations R{f; (p) in terms of ideal rotations R{r; cj)) as 



eral hundred neV. We also see that although the leading 
order of the error for the nine-piece pulse is higher than 
that of the other pulses, its coefficient is large enough 
that for f^ 2%, it is better to use the five- or 

seven-piece pulses. This suggests that generating similar 
pulses with even more pieces likely will not extend the 
range of values of dh for which low-error rotations are 
possible. 

For ho ^ Jmax, the noncommutation of the Hamilto- 
nian at different times even for Sh — makes the previous 
approach algebraically forbidding. However, in this case 
we have the possibility to perform rotations about axes 
f other than z. We make use of this freedom in order 
to construct SUPCODE pulses in a simple way: We first 
take an uncorrected (i.e., designed as if Sh — 0) rota- 
tion R{f; (f)) and then construct an uncorrected identity 



R{f; (j)) = R{f- 0) (/ + 5h A{f- 0) + 0{5hf) (3) 

and 

i{r;(l))^I + 6hB{f-^) + 0{6hf, (4) 

where A{f] (j)) and -B(f; (j)) are the first order corrections. 
For the corrected pulse M = R{l + Sh(A + B) + 0(6h)'^), 
it is clear that we want to make A{f; (f))+B{f; (p) — 0. Be- 
cause there are infinitely many ways to realize the iden- 
tity operation, there is sufficient freedom to engineer such 
a cancelation. For example, we consider a product of pa- 
rameterized forms of the identity such as 

ii{ai,bi) = R{hox + biJ^^^z; 2TTai) ■ R{x; 2tt) 

(5) 

• R{hax + biJ^^^^'z] 27r(l - a,)) 
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and 

l2{ai,bi,n) = R{x]2TTa.i) ■ R{hQX + hJ^nnxZ; 2n7r) 

(6) 

■Rix;2^^{l-a^)), 

where ai and bi are parameters satisfying < < 1, 
< 6i < 1 and n is a positive integer. 

As an explicit demonstration of SUPCODE, we show 
rotations around the x- and z-axes by arbitrary angles 
< (/) < TT for the particular case Jmax/^o — 4. (Com- 
binations of these are sufficient for universal one-qubit 
gates.) About x, the uncorrected rotation is performed 
by holding J = for a time (p/ho, and this is preceded 
by the identity / = /i(l/2, 61)2/2(1/2, 62, 15) where hi 
and 62 are chosen such that errors cancel (see Supple- 
mentary Information). About z, the uncorrected rota- 
tion is performed by a three-part puls d^^ l ^ ^ i?(z; (p) ^ 
R{x + z]Ti)R{x](j})R{x + z;7r), and this is preceded by 
the identity / = /i(l/2, &3)"^/2(a, ^4, 15) where a, 63, and 
&4 are chosen such that errors cancel (see Supplementary 
Information). As shown in Fig. [I] for = 7r/2, SUPCODE 
does indeed lead to a higher-order scaling of the error 
in (5/i/Jniax (and hence in fi/Jmax?2')5 ^'^d ^ reduction of 
error when (5/i/Jniax is less than a few percent. For in- 
stance, when (5ft./Jmax — 1%, errors are typically reduced 
by an order of magnitude. 

The tradeoff here is that SUPCODE rotations are typi- 
cally over an order of magnitude longer than uncorrected 
rotations. We note that for a given experimental set of 
parameters one should try optimizing the pulse sequence, 
which we have not done for the arbitrarily-chosen exam- 
ple above, as both the length and error of the corrected 
pulses could certainly be reduced by a significant con- 
stant factor by searching over different constructions of 



the identity, /. Nonetheless, for experiments in GaAs sys- 
tems with pumped nuclear spina^^ (/iq ~ 0.6/xeV, 5h ^ 
30neV) or micromagnet^ (/ig ~ 20neV, 5h ^ IneV) the 
sequences shown here could already deliver substantial 
improvement over the uncorrected ones. 

Realistic deviations from the ideal pulses assumed 
above would include charge noise, which adds a ran- 
dom quasi-static contribution to the exchange, and fi- 
nite rise times. The former could be reduced by using a 
multi-electron variant of the singlet-triplet qubilP^ and, 
in principle, it may even be possible to dynamically cor- 
rect by adding more degrees of freedom to the pulses. 
The latter could be compensated for by adjusting pulse 
parameters given the actual turn-on/off profiles of the 
pulse for a specific experimental setup. Thus, even in 
nonideal conditions SUPCODE could enable precise spin 
qubit rotations independent of shot-to-shot variation in 
the nuclear Overhauser field, also easing tasks such as 
ensemble-averaged measurements of singlet probability 
oscillations versus time by reducing hyperfine-induced 
decay. More importantly, this work allows satisfaction of 
the quantum error correction threshold within a substan- 
tially larger region of the physical parameter space than 
would otherwise be possible. The fact that T2 has now 
reached tens of microseconds in GaAs quantum dotJ^^, 
and millisecond^^ or even seconds in the presence of iso- 
topic purificatioiP^ in Si-based structures implies that 
gate errors are currently dominated not by dynamical 
fluctuations in the nuclear spin bath, but by the statisti- 
cal distribution of the magnetic field gradient, and these 
errors may be efficiently suppressed by SUPCODE. 
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Supplementary material 



In this supplementary material we provide detailed information related to the results shown in our main article 
"Composite pulses for robust universal control of singlet-triplet qubits". In Sec. |I]we give details of the expansion of 
the evolution operator with respect to fluctuations in the magnetic field gradient 5h, in the vicinity of — 0. Using 
these results we present, in Sec.|ll] a no-go theorem that no pulse sequence with a finite exchange interaction can be 
found to cancel all orders of error due to 5h. In Sec. |III| we show explicit forms of piecewise pulse sequences that, for 
/iQ = 0, cancels error due to 5h up to first, second and the third order in the evolution operator (which amounts to 
second, fourth and sixth order in the error per gate). In Sec. IV we show the pulse sequences that for ^ Q carry 



out X and z rotations while canceling the leading order error in Sh 



I. EXPANSION OF h AROUND ho = 



We start with the evolution operator defined in the main text 







|exp 









with 



J{t) 



h — ho + 6h. 

Expanding the evolution operator in powers of Sh in the vicinity of ho = 0, 



(S-1) 



(S-2) 



(S-3) 



n=0 



one can show that 



n 




fir) 



fc=i 



An + sin 



fir) 



fe=i 



Bn 



(note that the product is in descending order in observance of the time-ordering of operators) where 



f{T) = / dtJ{t), 
Jo 



(S-4) 
(S-5) 

(S-6) 



A-n = 



n even 



CTj; n odd 



Bn — 



—i(Tz n even 
ay n odd 



(S-7) 
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and / is the 2x2 identity matrix and a^, cry, and cr^ are Pauli matrices. 
Thus, for the nth-order term to vanish, one must have 



n 



t' 



dtL 



exp 



(s- 



k=l 



II. NO-GO THEOREM FOR CANCELING ERROR TO ALL ORDERS 



In this section we show that there does not exist a pulse sequence with a bounded exchange interaction that 
suppresses all orders of error due to 6h. 
Consider a function 



V'(T,A)=i+5]A" n / 

71 — 1 \7n— n ^ 



dt'ra exp 



I 

L k=l 



+ \ [ dt[ exp [-ifit[)] +\^ [ dt'^ [ ' dt[ exp + 
Jo Jo Jo 

+ A3 r dt', dt', /*' dt[ exp l-ifiQ + if{Q - ifiQ] + ■■■ 



JO JO 

which solves the partial differential equation 



with boundary condition 



dip 
df 



= Xexp[-tfiT)]r, 



V'(r = 0,A) = l. 



(S-9) 

(S-10) 
(S-11) 

(S-12) 
(S-13) 



Observe that if we can find a pulse which cancels all orders in Eq. (S-8) at T = Ty, it means that Eq. (S-12 1 has a 
solution that also satisfies 'ip{T — Tf, X) ~ 1 for any A. To simplify the problem further, we write 



7/;(r, A) = exp 



:f{T) 



exp 



u{T, X) — i 



w{T,X) 



where u{T, A) and w(T, A) are real functions satisfying 

du{T, A) 



dT 
dw{T, A) 
df 



Acosw(T, A), 
-2Asinw(r, A) - J{T). 



The associated boundary conditions are 

u(r = 0,A) =0, 



u{T^Tf,X)^X dT' cos w{T',X) = 0, 
Jo 

w{T = 0,X) =0, 
wiT^Tf,X) = -f{Tf). 



(S-14) 

(S-15) 
(S-16) 

(S-17) 
(S-18) 

(S-19) 
(S-20) 



exists in this limit and is given by 



We focus on Eqs. (S-16), (S-18l-(S-20) only and consider the limit A — >• cx). It is clear from Eq. (S-16) that w{T,X) 



(S-21) 



whe re n i s an integer. To satisfy the boundary condition in Eq. ( |S-19 ), we need n = and J(0) = 0. To satisfy 
Eq. (S-20), we must have 



w{Tj,X)^~^^-f{Tj), 



(S-22) 
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which in turn requires J{Tf) = f{Tf) = 0. This already negates the possibihty of performing any nontrivial {f{Tf) ^ 
0) z-rotation. However even error- free trivial rotations are impossible when the final condition, Eq. (S-18), is taken 
into account since this immediately gives T/ — 0. We conclude that it is impossible to design a J(t) that cancels 
field gradient errors to all orders. In reaching this conclusion, we have implicitly assumed that J(t) is bounded and 
T/ < cx). 



III. PIECEWISE PULSE SEQUENCES FOR /lo = 

The no-go theorem shown in Sec. |ll] implies that it is instructive to design pulses which cancel successive orders. In 
this work we focus on the piecewise constant pulse, which can be expressed as 



' 6h Jfc ^ Tk 



2 2 / J, 



max 



(S-23) 



k=N 

where we have defined the dimensionless quantity — Jmax^fe for convenience. 

Let us also define some simplifying notations, tk refers to the duration of a pulse on the k*^ piece. The total 
duration of time after the fc*'' piece would be 

fe 

Tk^Y.^,, (S-24) 
with To = and Tjv = Tf indicating the initial and final time. J{t) can be expressed as 

N 

J{t) = J2 '^'^^(^ - Tk-i)Q{Tk - t), (S-25) 

k=l 

where Q{t) is the Heaviside step function. 

The motivation is to cancel successive orders of 6h in the expansion with number of pieces N as small as possible. 
We will impose further constraints to simplify the search: First, we only allow < Jfc < Jmax- The non-negativity 
of this prescription ensures that the designed pulse is relevant to experiments on singlet-triplet qubits. Second, we 
consider symmetric pulses, i.e., J{t) = J{Tf — t), in this section. This ensures that U{Tf,0) has no ay componenlP. 



A. Three-piece pulse sequence 



One of the simplest ways to cancel the leading order error with strictly positive values of J is via a three-piece pulse 
sequence. In Eq. (S-23) we take = 3 and the three-piece pulse sequence can be characterized by 



(ti, Ji) = (rs, J3) = (00, Jmax); (t2, J2) = (47r - 2(f)Q, Jmax/2), 



(S-26) 



where (f)o is the desired rotation angle around the z-axis. It is straightforward to verify that the evolution operator 
under this pulse is 



UiTf,0) 



exp 



■ 90 

2 



Sh^ 
2 p 



(47r - 00 + sin 00 ) 



sin — / 
2 



I cos (T, 

2 



^ J max ^ 



(S-27) 



Here the deviation from the desired rotation has been suppressed up to second order in 8h. Note that this pulse 
actually sweeps the Bloch vector through an angle = 27r -I- 0Oi which is the origin of the trivial additional phase 
factor. It is clear that in order to achieve error cancellation, the Bloch vector generally must be swept through more 
than 271" about the z-axis since its path is deflected from the ideal {8h ~ 0) path in opposite directions in the "eastern 
and western hemispheres." [One can also explicitly show the necessity of larger angles from Eq. (S-4).] Thus the 



three-piece pulse sequence cancels leading order error simply by ensuring that, in the absence of ^/i, the Bloch vector 
spends an equal amount of time in each hemisphere during its rotation. In this section, we will use 0o to denote the 
"desired" rotation angle, and the "actual" rotation angle, where is larger than 0o by a multiple of 27r. 
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B. Five-piece pulse sequence 



In this subsection we set = 5 in Eq. ( |S-23 1, which allows us cancel the dependence of U{Tf,0) on 5h up to the 
second order, using a symmetric pulse with parameters ri = T5, T2 = T4 = 0/4, = J3 = J5 = 0, J2 = J4 = Jmax- 
We expand C/(T/,0) as 



C/(r/,0)= ( 



= I cos-/ - i sin -cr. 



Sh 



/ (f> (f>\ 

+ 2ti cos - + 2 sm - J cr^j 



8 P 

max 



4tiT3 + (4ti^ + T^) cos - + 2(4ri + 2t3 + 0) sin [4t3 - 2(2r3 + (f)) cos 2 + (4 + ^"3 ) sin ^ 



(/ If 



(S-28) 



To make the first order coefficient vanish, one must choose 



Ta -2 ( Ti cos 2 + ^™ 2 



(S-29) 



Plugging Eq. (S-29) into Eq. (S-28 1, it suffices to satisfy 



-2ri - + 4ti cos - - 2ri cos (/) + 4 sin - + {t'^ - I) sin (/> = 



(S-30) 



to make the second order terms vanish. We choose a root that is positive for 2tt < cj) < Stt, that is 



Ti = CSC (/) I I — 2 cos — -I- cos (j) + y 4: ~ 8 cos 2 ^ '^'^^ ' 



(S-31) 



To summarize, to achieve a R{z;(j)) rotation while canceling the dependence of U{Tf,0) on Sh up to the second 
order, we should choose the pulse parameters as 



T5 = Ti — CSC ( I — 2 cos — -I- cos ( 



T4=T2= (/>/2, 

T3 = -2 ( Ti COS - + sm ^ 



4 — 8 COS + 4 cos (j) + (j> sin < 



(S-32) 

(S-33) 
(S-34) 



A plot of these parameters as a function of cf) is given in Fig SII Note that this pulse is defined for 4> & (27r, 37r), 
which is equivalent to a rotation R{z; 0o) with (po g (0, tt). Rotation of angles outside the range (0, tt) may be achieved 
by duplicating existing pulse sequences. For example, to achieve a £-rotation of (f>Q = 1.27r one could apply twice the 
00 = O.Gtt (corresponding to = 2.6n) rotation. 

As — >■ 2iT, Ti — >■ 00, T2 = 7r/2, Ta 2ti. If we want to fix the total duration of the sequence, we let Jmax — >■ 00. 
The pulse sequence becomes the CPMG puls^. This is no surprise. In this case, it can indeed be shown that for a 
pair of instantaneous tt pulses, setting leading order errors to zero while maintaining the time-reversal symmetry of 
the pulse sequence enforces the Uhrig conditioiP. 



C. Seven-piece pulse sequence 



In Eq. (S-23) we take N = 7 and ti — T7, T2 = T4 = rg = 0/3, T3 = rs, Ji = J3 = J5 = Jj = 0, J2 = Ja = Je = Jmax- 
Since the number of independent variables are the same as for the five-piece pulse, it can still cancel the dependence 
of U{Tf,0) on Sh up to second order, but not to third order. However, it will expand the range of rotation angles for 
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1 — ri — r-^T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




Q I I I I I I I I I I I I I I I I I I I i_ 

271 2.2k 2.47t 2.67t 2.871 371 

(t)[rad] 



FIG. SIl: Parameters for the five-piece pulse as defined in Sec. |III B] 



which the pulse can directly be applied. We expand U [Tj , 0) as 

6h 



U[Tf,0) = l^cos -I -ism -(7,_ j - I J [T3COS - + Ti cos - + sm -jfi. 



4 /2 



r3(4ri + T3) cos ^ + (2rf + Tg^) cos ^ + (0 + 4ti + 4T3) sin | 



4t3 cos ^ - (4r3 + (/)) cos ^ + r| sin ^ + (t| + 2) sin ^ 



f 5h \3 



Similar to the previous section, to make the first order coefficient vanish, one must choose 



T3 = - sec - l^Ti cos - + sm - j . 



Plugging Eq. (S-36) into Eq. (S-351, it suffices to satisfy 



20 (f) 
8ti + d + 4ri cos — 3 cos — 



d) (j) . n . 5(p o 7(j> 

i COS ^ — 6 cos — — (tj — 5) cos + ('''1 ~ 1) cos — 6 sec 



to make the second order terms vanish. We choose a root that is positive for 3tt < (f) < Stt, that is 



(S-35) 



(S-36) 



sin — = 
2 



2 cos t — cos t — cos ^ + 2 

2 6 6 



COS ^ — 6 sin ^ j f sin | — 2 sin ^ 



sm 7? — sm -if 

6 6 



(S-37) 



(S-38) 



To summarize, to achieve a R{z; <f) rotation while canceling the dependence of U(Tf, 0) on Sh up to second order, 
we should choose the pulse parameters as 



T7 = n = 

T4=T2= (jy/i, 



2cos I - cos I - cos ^ + 2Wcos2 I Ucos f - 6sin f 1 (sin f - 2sin g 



• (b ■ 5<i 

sm ^ — sm 

6 6 



(ncos-+sm-J. 



Ts = r3 = - sec - yri cos - -1- sm ^ 



(S-39) 

(S-40) 
(S-41) 
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T — n — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



Xt 




Q I I I I I I I I I I I I I I I I I I I l_ 

371 3.4JI 3.8JI 4.271 4.671 571 

(t)[rad] 



FIG. SI2: Parameters for the seven-piece pulse as defined in Sec. |III C[ 



A plot of these parameters as a function of (j) is given in Fig SI2 Note that this pulse is defined for <f> € (Stt, Stt), 
which is equivalent to a rotation R{z; (po) with (f)Q S (— 7r,7r). 

As (/) — >■ Stt, both Ti, T3 — )• oo, T2 — — TT. However, in this limit, we have 



Ii^^i±^.l + V2. 
Ti 2 + ^/2 



(S-42) 



Note that sin | = ^2 - ^2/2, and 



sin^ ^ - sin^ f ^ 

^ = 1 + ^2. 

sin^ I 



(S-43) 



Therefore in this limit the pulse is nothing but the n = 3 Uhrig pulse, which can be understood in a similar way as 
in the previous section. 



D. Nine-piece pulse sequence 



In Eq. (S-23) we take iV = 9 and the nine-piece pulse sequence can be expressed as ti = rg, T2 = T4 = tq = — 0/4, 
T3 = T7, Ji = J3 = J5 = J7 = Jg = 0, J2 = Ja = J& = Jf, = Jma.x- Now there is one more independent variable, 
which means that we are able to cancel the dependence of U{Tf,0) on Sh up to the third order. The expansion of the 
evolution operator is complicated, so we start with the first order in Sh, which is 



—I 



Sh 



+ 2t3 cos — -I- 2ti cos ^ + 2 sin — J a^- 



(S-44) 



Making this term vanish requires 



T5 = -2 ( T3 cos ^ + Ti cos - + sin - 



(S-45) 



Using Eq. ( S-45 ) , the expansion simplifies to a function only dependent on ri , T3 and 4>. To further simplify the 
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FIG. SIS; (Color online.) Parameters for the nine-piece pulse as defined in Sec. HID The area shaded by light yellow indicate 
the regime where positive ri cannot be found. 



notations we introduce functions 



/i(n,T3, 



/2(ti,T3,, 



3<p 



— 2ri — 4t3 — (p + 6r3 cos — + 4ti cos — — 2t3 cos ^ 



3(/) 



2ti cos (j) - 2riT3 sin ^ + (rf + 4) sin - + 2riT3 sin + (t^ - 1) sin 0, 



2t? cos ■ 



+ 3ri(3 + rf) + 3(4t3 + cos | - 2r3(15 - 3rf + t|) cos ^ 



6ti (t| + 4) cos (j) - 6{t^ - 1)t3 cos ^ - 2ti {t^ - 3) cos y - 3r3(4T3 + (/)) sin ^ 
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- 1 + - 3t| + Ti (4r3 + 0) sin ^ + 36riT3 sin y + 6{2t^ - t| - 2) sin </> 



- 12tiT3 sin y + 2(1 - 3rf ) sin y . 
With these definitions, we express the expansion of the evolution operator as 



(S-46) 



(S-47) 



U{Tf,0) = [ cos J/-isin Ja, 



4 P 



o 



-/i(n,T3,( 
Sh 



sin — / + i cos —a^ 



' max 



^24J3 



( — ) 

^ max ^ 



(S-48) 



To cancel the dependence of U{Tf, 0) on 5h up to third order, we must find ri and as positive, real solutions to 
the coupled nonlinear equations 



/i(ti,T3,0) =0 

/2(Tl,r3,0)=O' 



(S-49) 



while at the same time requiring a positive T5. 

We found a numerical solution to the above equations for (j) S (47r, 4.56547r)n(57r, Gtt), which corresponds to rotation 



R{z;(j)a) with 0o € (0,0.56547r) H (tt, 27r). The numerical solution is shown in Fig. SI3 



IV. PIECEWISE PULSE SEQUENCES FOR ho ~ J^^^/A, 

As described in the main text, we perform corrected rotations as RI where the implementation of the identity / is 
chosen such that the first order errors in 5h exactly cancel those from the rotation R. Our implementations of / are 
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h h R 




100 150 

ax 




(b) 



0.2re OA-n 0.6Jt 0.8Jt 
^ [rad] 



FIG. SI4: (Color online.) Panel (a): Example of a pulse sequence for rotation about the a;-axis by 7r/2 with /iq = Jmax/4. 
Panel (b): Parameters for rotations about the a;-axis for the same ho value and < </> < tt. 




(a) 

FIG. SI5: (Color online.) Panel (a): Example of a pulse sequence for rotation about the z-axis by 7r/2 with ho = Jmax/4. Panel 
(b): Parameters for rotations about the ^-axis for the same ho value and Q < cj) < n . 



of the form 



(S-50) 



where the parameters Ci, hi are restricted by the physical constraint < J{t) < Jmax- For the specific forms Ii, I2 
given in the main text, this translates to the constraints 



< a, < 1, 
< h < 1. 



(S-51) 
(S-52) 



The particular implementation of / that will be chosen for a given experimental situation will depend on the par- 
ticular rotation R as well as the experimental parameters ho/ Jmax, with a combinatorial search necessary to find an 



implementation that: (a) cancels the first-order error from R; (b) satisfies the constraints of Eqs. (S-51) and (S-52 1; 
and (c) is close to 'optimal' in some experimentally meaningful sense, such as having the overall shortest duration 
in time, having the smallest second-order error in Sh, or being least sensitive to over-/under-shoot in J{t). As an 
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example of our technique we found sequences valid for ho ~ Jinax/4, for corrected rotations about angles < < tt 
around the x- and z-axes. For the x-axis rotation we have 

R{x; (j)) + 0{Shf = R{x; 0)/i(l/2, 61)2/2(1/2, 62, 15). (S-53) 

The first-order term in Sh is zero when bi^2 are chosen to satisfy 



- 3, (S-54) 

2^h^ + bl + h ^ 15^2/^ , . 

{h^ + blf 2b,{h^ + bl)V^' ^" ^ 

where we have introduced the dimensionless parameter h = /io/>Anax- The solution to these equations for /i = 1/4 are 
shown in Fig. SI4 In Fig. |SI4[ a) we show a specific example of the pulse sequence for 7r/2 rotation about the a:-axis. 

Fig. |SI4[ b) shows solutions for a range of rotation angle < (/) < tt. For certain values of close to tt, bi exceeds 1 
by up to 7% but this is not a problem because Jmax is not typically a hard limit for pulse amplitudes. 
For the same parameter range we construct the corrected z-axis rotation as 

R{z; (j)) + 0{5hf = R{x + z; n)R{x] (j))R{x + z; 7r)/i(l/2, bifi2{a, 64, 15), (S-56) 

with a and 63^4 satisfying 

Qh^ iGh? Up siufj) cos(/) + l_ 

{h^ + bif ^ {h^ + bif ^ + ^ + - ' ^'"''^ 

cos = 30^^4/.^sin(2^a) ^ ™0 ^ ^ (S_58) 
{P + blf 2^/2 
3064/icos(27ra) I263 663/1 <^ 1 n 



(^2+52)3/2 ^^2+52 (^2+^,2)3/2 ^J^ 



(S-59) 

The solution to these equations is shown in Fig. |SI5[ We show in Fig. |SI5[ a) an example of the pulse sequence 
appropriate for rotating about the z-axis by 7r/2, while in Fig.ISI5[b) we show solutions hi < (j)< n. 



^ To see this fact, note that the operator exp[~i{hax/2 + Jkaz/2)tk] can be written in the form Ai — Uiol + ai^a^ + ai^a^ with 
a-io, o.ix, diz arbitrary complex numbers. It is then straightforward to show that for any operators Ai and A2 with arbitrary 
coefficients, A2- A\- A2 can also be written in such a form, free of ay terms. Applying this statement recursively to Eq. ( S-23 1 , 
one sees that for any J{t) satisfying J{t) = J{Tf — t) the resulting evolution operator U does not contain ay component. 

^ G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 

^ H. Y. Carr and E. M. Purcell, Phys. Rev. 94, 630 (1954); S. Meiboom and D. Gill, Rev. Sci. lustrum. 29, 688 (1958). 



